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Abstract 

Proposals for quantum computing using rotational states of polar molecules as qubits have previously 
considered only diatomic molecules. For these the Stark effect is second-order, so a sizable external electric 
field is required to produce the requisite dipole moments in the laboratory frame. Here we consider use 
of polar symmetric top molecules. These offer advantages resulting from a first-order Stark effect, which 
renders the effective dipole moments nearly independent of the field strength. That permits use of much 
lower external field strengths for addressing sites. Moreover, for a particular choice of qubits, the electric 
dipole interactions become isomorphous with NMR systems for which many techniques enhancing logic 
gate operations have been developed. Also inviting is the wider chemical scope, since many symmetric top 
organic molecules provide options for auxiliary storage qubits in spin and hyperfine structure or in internal 
rotation states. 
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I. INTRODUCTION 



In principle, a quantum computer can perform a variety of calculations with exponentially 



fewer steps than a classical computer 



jg]. This prospect has fostered many proposals for means 



to implement a quantum computer II Vl-4 1 VTI . Using arrays of trapped ultracold polar molecules is 



considered a promising approach, particu 



obtain large networks of coupled qubits 1115 



ar 



y since it appears feasible to scale up such systems to 



2911 . Molecules offer a variety of long-lived internal 



states, often including spin or hyperfine structure as well as rotational states. The dipole moments 
available for polar molecules provide a ready means to address and manipulate qubits encoded in 
rotational states via interaction with external electric fields as well as photons. 

Entanglement of qubit states, a major ingredient in quantum computation algorithms, occurs 
in polar molecule arrays by dipole-dipole interactions. In a previous study, we examined how the 
external electric field, integral to current designs for quantum computation with polar molecules, 
affects both the qubit states and the dipole-dipole interaction [|29j]. As in other work concerned 
with entanglement of electric dipoles, we considered diatomic or linear molecules, for which the 
Stark effect is ordinarily second-order. Consequently, a sizable external field (~several kV/cm) is 
required to obtain the requisite effective dipole moments in the laboratory frame. 

In considering the operation of a key quantum logic gate (CNOT), we evaluated a crucial pa- 
rameter, Ao), due to the dipole-dipole interaction. This is the shift in the frequency for transition 
between the target qubit states when the control qubit state is changed. For candidate diatomic 
molecules, under anticipated conditions for proposed designs, Ao) is very small (~20-60 kHz). 
It is essential to be able to resolve the Ao) shift unambiguously, but in view of line broadening 
expected with a sizeable external field, whether that will be feasible remains an open question 



This question led us to consider polar symmetric top molecules, for which the Stark effect is 
first-order in most rotational states. The effective dipole moments are then nearly independent 
of the field strength. That enables use of a much lower external field (a few V/cm) to address 
and manipulate the dipoles, improving prospects for resolving the Ao) shift. The constancy of 
the symmetric top effective dipole moments also makes entanglement properties of electric dipole 
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interactions isomorphous with those for nuclear magnetic resonance systems. This suggests that 
NMR techniques, extensively developed for quantum computation but limited in application by the 



small size of nuclear spins and scalability prospects 



17. 



30 



3 lh might find congenial applications 



with qubit systems comprised of polar symmetric top molecules. 

H. EIGENSTATES FOR A POLAR SYMMETRIC TOP 

The Hamiltonian for a single trapped polar symmetric top molecule in an external electric field 
may be written 

H = H R + H s +H r + H,, ? (1) 
The major term is the rotational energy 

H R = BJ 2 + (A - B)J 2 Z (2) 

where J denotes the total rotational angular momentum and J z its projection on the symmetry axis; 
A and B, the rotational constants, nominally inversely proportional to the moments of inertia about 
the principal axes along and perpendicular to the symmetry axis, respectively (actually effective 
values averaged over vibration and centrifugal distortion of the molecule). The Stark energy from 
interaction with the external electric field is 

H 5 = -fx • £ = -//8cos# (3) 

with 6 the angle between the body-fixed dipole moment /i (along the symmetry axis) and the 
direction of the field. The trapping energy is 

P 2 

Ht = 2^ + V,mp (4) 

but at ultracold temperatures the translational kinetic energy p 2 /2m is quite small and very nearly 
harmonic within the trapping potential V, rap ; thus H T is nearly constant and for our purposes can 
be omitted. The remaining term, H M , represents interactions arising from nuclear spins and/or 
quadrupole moments; here we omit treating these, except for an important effect of the quadrupole 
interaction in modifying selection rules for transitions between qubit states. 



In familiar notation, 021 13311 the eigenenergy for H R is 



E R (J, K) = 57(7 + 1) + (A - B)K 2 



(5) 



For a prolate top, A > B; for an oblate top, A < B. The Stark energy for H 5 is 



E s (J,K,Mj) = -fi£MjK/J(J + 1) 



(6) 



to first order. The second-order term is far smaller B4I1 (so neglected here) except for K = or 
Mj = states (which we will not use as qubits). The corresponding eigenfunction for H s can be 



written as 11321 . 



33] 



\J,K,M) = {-\) 



M-K 



27+ 1 
~8n 2 ~ 



e m d ] _ M _ K (ey* K 
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where (p, 9 and x are the Euler angles and d J M K {6) is a Jacobi polynomial (aside from a simple 
prefactor). Hence, in addition to the polar angle 6 that governs the Stark interaction, the eigen- 
function depends on the azimuthal angles^- and <p associated with, respectively, the projections of 
J on the molecular symmetry axis and on the B-field direction. 

Figure CD displays for the K = 1 sublevels of the 7=1 and 2 symmetric top rotational states 
the (a) eigenenergies W = E R + E s and (b) expectation values (cos) = //e/y/yu for the projection of 
the dipole moment on the field direction, as functions of fiE/B. The dependence on fiE/B differs 
markedly from a similar plot for a diatomic molecule (for which K = 0; cf. Fig. 1 of ref. 112911 ): 
there the effective dipole moments are field-dependent and vanish at zero-field. For symmetric top 
qubit states, to take advantage of the first-order Stark effect, we consider only K ± and Mj ± 
states. For such states, the effective dipole moments, 



fi eff = -dEs/dE = /uMjK/J(J + 1) 



(8) 



are just constants independent of the field (except at unusually high fields, where higher order 
terms become important [34]). According as /j. e ff is positive or negative, the Stark energy drops or 
climbs as the field strength grows, so the molecular states are termed high field seeking (HFS) or 
low field seeking LFS), respectively. 
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A. Choice of qubit states 



We consider two qualitatively distinct choices for qubit states, designated I and II. The orthodox 
choice, type I, is exemplified by 

|0> = |/= l,K= l,Mj = -1> and |1> = \J = 2,K = l,Mj = -1) (9) 

For this choice (green curves in Fig. Q]), radiation induced transitions between the qubits are fully 



allowed, in accord with the familiar selection rules, A J = 0, +1; AK = 0; AM = 0, +1 113211 . Also, 
both the |0) and |1) qubit states are LFS, thereby facilitating trapping by either DC or AC fields or 
an optical lattice [35]. The corresponding eigenenergies, E R + E s , are 

W =A + B + ^f- and W l = A + 5B + ^ (10) 
2 6 

and the cos# matrix elements are 

Co = <0|cos#|0> = -\, d = <l|cos0|l> = ~, C x = <0|cos#|l> = (11) 

I o ID 

We are particularly interested in an unorthodox choice, type II (red curves in Fig. [[])• For this, 
the qubit states are 

|0> = \J = l,K = l,Mj = +1) and |1> = \J = 1,K = l,Mj = -1> (12) 

The eigenenergies are degenerate at zero-field but for £ > split apart strongly and linearly, 

W Q =A + B-^ and W X =A + B + ^- (13) 

and the cos# matrix elements are 

Co = <0|cos#|0> = 1 d = <l|cos0|l> = -I C x = <O|cos0|l) = O (14) 

These type II qubits render the effective dipole moments constant and equal in magnitude but 
opposite in sign. However, type II qubits require further specification. As initially defined in 
Eq.(12), the transition |0) <-» |1) between the qubits requires AMj = ±2. Thus, it is not allowed 
as a one-photon electric dipole transition (the transition cosine, C x = 0). It is allowed as a two- 
photon transition (using the J = 1, K = 1, Mj = state as intermediate). Another remedy, 



simpler to implement, is to use a molecule that contains a nuclear quadrupole moment. Even a 
small quadrupole coupling constant typically introduces sufficient mixing of Stark states to make 



AM j = +2 transitions become prominent in microwave or radiofrequency spectra 113 611 . In accord 



with theory Q32 



3711 . in the next subsection we show that modifying the type II qubit choice to 
exploit the quadrupole hyperfine structure renders C x + 0, enabling |0) <-» |1> to occur as a one- 
photon transition. 

In another contrast with type I, for type II qubits |0) is HFS while |0) is LFS. That is also 
often the case for qubit states considered for diatomic molecules, and is not regarded as a seri- 



ous handicap fl33fl - Although HFS states are harder to trap, both HFS and LFS can be captured 



simultaneously in an AC trap or an optical lattice 



B. Quadrupole perturbation of Stark states 



For simplicity, we consider symmetric top molecules having only one atom with a nuclear 
quadrupole moment, with that atom located on the symmetry axis. We also treat explicitly only 
cases in which the nuclear spin 1 = 1 for that atom, and the quadrupole interaction is much smaller 



than the Stark energy. The CH^CN molecule [|38|l is a prototypical case: for the 14 /V nucleus (spin 



1=1), the quadrupole coupling constant is eqQ = -4.22 MHz. For conditions in prospect for a 



quantum computer, usually //£ > 



00 MHz. A first-order perturbation treatment, referred to as the 



3711 . and governed by the ratio eqQln&, hence is appropriate for 



"strong-field approximation" 11321 
this example and many others. 

When set-up in the usual \JKMjIMj) basis, with Mj the projection of the nuclear spin on the 
E-field direction, the Hamiltonian matrix, H R + H s + H G , is diagonal in J, K, and /. The H R and 
H s portions are also diagonal in Mj and Mj whereas H e has off-diagonal elements which connect 
Mj and Mj states differing by up to two units. In consequence of the resulting mixing, neither Mj 
nor Mj is a "good" quantum number. Their sum, Mj + Mj remains good, however, since the total 
angular momentum along the field must be constant. Accordingly, we modify our choices for the 
|0) and |1) qubits of Eqs.(9) and (12), that involve Mj = +1, to specify them further as particular 
hyperfine components with Mj + M/ = 0. In Appendix A we evaluate the contributions from H e 
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to the qubit eigenenergies and cosine matrix elements. 

In first-order, the quadrupole interaction simply adds to the qubit eigenvalues of Eq. (10) or 
(13) a diagonal term given by 

Eq = eqQ/40 or eqQ/56 (15) 

for J = 1 or J = 2, respectively. 

The cosine matrix elements of Eqs.(ll) and (14) are augmented by terms involving w = 
\eqQ\/li£, given in Table I. Since typically w < 10 2 , these contributions are insignificant for 
type I qubits, and for the Co or C\ elements for type II qubits, but of major importance in the Cx 
transition element for type II, which would otherwise be zero. Even when Cx is very small, con- 
ventional power levels suffice to make transitions facile between the Mj = +1 Stark components 

y. 

TABLE I: Cosine matrix elements for symmetric top qubits". 

Type I qubits Type II qubits 

Co -l/2-0.00168w + 0.0418w 2 1/2 - 0.00347 w - 0.021 3w 2 

Ci -1/6 + 0.00526w + 0.02 18w 2 -1/2 - 0.00168w + 0.0418w 2 

C x Vl5/10-0.00658w -0.0437w 2 + 0.153w - 0.0108w 2 

"Terms in w = \eqQ\//j£, are contributions from quadrupole coupling. These were fitted to 
results of numerical calculations (see Appendix A) extending over the range w < 1. 



III. TWO INTERACTING DIPOLES 

Adding a second trapped polar symmetric top, identical to the first but a distance r l2 apart, 
introduces the dipole-dipole coupling interaction, 

Hi fi 2 -3(Mi -n)0w 2 -n) 

Vd-d = : 1, (16) 

|ri - r 2 | 3 

Here n denotes a unit vector along r 12 . In the presence of an external field, it becomes appropriate 



to express V^-d m terms of angles related to the field direction (Appendix A in ref. I129h ). The 
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result after averaging over azimuthal angles reduces to 



Vd-d = £2(1 - 3cos oOcos^cos^ 



(17) 



where £2 = fi /r n , the angle a is between the vector and the field direction and polar angles 9\ 
and 8 2 are between the /ii and fi 2 dipoles and the field direction. 

When set up in a basis of the qubit states (either type I or II) for the pair of molecules, 
{|00>, |01), |10), |11)}, the H R + H s portion of the Hamiltonian takes the form 



Wq + W q 

o w + w: 



o 












Wi + W' Q 

Wi + W[ 



and the V d ^ d portion is 



(18) 



(19) 



CoCq CqC' x CxC' q CxC' x 
CqC' x CoCj CxC' x CxC\ 
CxC'q CxC' x C\C' q C\C' x 
<CxC' x CxC[ C\C' X C\C[ t 
where £2„ = £2(1 - 3cos 2 Q').The primes attached to quantities for the second dipole indicate that 
the external field magnitude will differ at its site; that is necessary for addressing the sites and to 
ensure that the qubit states |01) and |10) differ in energy. 



A. Evaluating entanglement of eigenstates 

The form of the Hamiltonian in Eqs. (18) and (19) is identical to that for two polar diatomic 
molecules, treated in ref. Thus, we follow the same procedures in evaluating eigenstate prop- 
erties and entanglement for symmetric tops, merely introducing the appropriate matrix elements 
for qubits of types I and II (as specified in Sec IIA). We again use unitless reduced variables, x = 
fiE/B and y = £2 a /5; in terms of customary units, these are given by 

x = fiE/B = 504/z(Debye)8(kV/cm)/ J B(MHz) (20) 
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y = QJB = 1.51 x KrV(Debye)/r 3 (//m)/fl(MHz) (21) 

Likewise, we use z = eqQ/B for quadrupole coupling terms. The pertinent ranges are x < 1, y < 
10~ 5 , and |z| < 5 x 10~ 3 for candidate symmetric tops (with dipole moments n < 4 D, quadrupole 
coupling \eqQ\ < 10 MHz, and rotational constants B > 2000 MHz) under conditions deemed 
practical for prospective quantum computer designs (field strengths £ < 1 kV/cm, intermolecular 
spacings r ~ 0.5/im). Unless otherwise noted, we take a = 90°. In the pertinent regime, the 
dependence on x, y, and z of the eigenenergies is simply linear in all three variables. 

Another key variable is Ax = x' - x, specified by the difference in the field strength at adjacent 
qubit sites. As the site addresses are provided by observing the one-qubit transition, |0) <-> |1), the 
size of Ax must be large enough to produce a clearly resolvable Stark shift between the sites. Yet 
Ax must not exceed X R /N, where N is the number of sites and X R the range in x of field strengths 
considered feasible. To benefit from keeping the field strength relatively low, we take X R ~ 1 ; 
then to accommodate N sites requires Ax < X R /N. At least for exploratory calculations for up to 
N ~ 10 3 , we consider 10" 4 < Ax < 10" 2 appropriate. 

Tables II and III exhibit properties, for qubit types I and II, respectively, of the four eigenstates 
of the two-dipole system, listed in order of increasing energy (z = 1 — » 4). The eigenvalues are 
obtained as simple explicit functions of x, x' , y, z, applicable to any polar symmetric top molecule 
and conditions within the pertinent regime specified above. Also indicated, in order of magnitude 
only, are quantities that express the extent of entanglement among the qubit basis states, but must 
be evaluated by numerical means. Entanglement is exhibited most directly in the coefficients with 
which the qubit basis states appear in the eigenfunctions, 

= fljIOO) + 6,|01> + c/110) + di\U) (22) 

In Appendix B, we give somewhat cumbersome formulas for these coefficients in terms of x, x', 
y, z. Tables II and III show just orders of magnitude, evaluated for CH 3 CN, under conditions 
specified in Table IV. This is done to illustrate most simply a major point: In the pertinent range, 
the entanglement is so feeble that the successive eigenfunctions *P; differ only slightly from the 
respective basis qubits, {|00>, |01), |10), |11)}; there is little admixture with other qubits. 
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TABLE II: Eigenproperties for N = 2 symmetric top dipoles, type I qubits a . 



(Ei-2A)/B |00> |01 > 1 10) |11> C 



12 



1 2 + ! + T + 4- + 2% 1-O(10~ 15 ) +<9(10- 8 ) +<9(10^ 8 ) -6>(10^ 8 ) <9(10- 8 ) 

2 6 + f + f + ^ + |§ -<9(10~ 8 ) 1-O(10^ 7 ) -0.0009 +<9(10- 8 ) 0.0018 

3 6 + ^ + f + ^ + -<9(10" 8 ) -0.0009 1 - <9(10- 7 ) +<9(10^ 8 ) 0.0018 

4 10 + f + £ + ^ + ^ +O(10" 8 ) -<9(10~ 8 ) -<9(10^ 8 ) 1 - O(\0~ 15 ) <9(10~ 8 ) 
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a Herex = = 0.0107, y = Q a /B = 2 x 10~ 6 , z = eqQ/B = 5 x 10~ 4 , Ax = x' - x = 10 



TABLE III: Eigenproperties for N = 2 symmetric top dipoles, type II qubits". 



i 


(E r 


- 2A)/S 




100) 


101) 


|10) 


111) 


Cu 


1 


2- - - 

* 2 




20 


1 - O(10" 17 ) 


-<9(10~ 9 ) 


-<9(10^ 9 ) 


-<9(10~ 12 ) 


O(10" 12 ) 


2 


2 + - 


x' y 

2 4 ^ 


z 

20 


+<9(10~ 9 ) 


1 - O(10" 17 ) 


-<9(10^ 10 ) 


+<9(10~ 9 ) 


<9(10~ 9 ) 


3 


2+ * 


x y 

2 4 ^ 


z 

20 


+<9(10" 9 ) 


+<9(10~ 10 ) 


1 - COO" 17 ) 


+<9(10" 9 ) 


<9(10" 9 ) 


4 


2 + f + 


y y 


z 

20 


+<9(10~ 12 ) 


-<9(10~ 9 ) 


-<9(10~ 9 ) 


1 - <9(10^ 17 ) 


0(1O~ 12 ) 



fl Here x = fi£/B = 0.0107, y = Q a /B = 2 x 10~ 6 , z = eqQ/B = 5 x 10 -4 , Ai = i'-i= 10" 



TABLE IV: Parameters for CH 3 CN molecule. 



Properties 
fx 3.92 D 

B 9198.8 MHz 
eqQ -4.22 MHz 
yu£ 988 MHz 
Q a 18.5 kHz 



Reduced variables" 
x = fxE/B = 0.0107 
Ax = //(8'-8)/fl - 10" 3 
y = Q a /B = 2 x 10~ 6 
Z = eqQ/B = 4.6 X 10~ 4 
w = \eqQ\/n£ = 4.3 X 10~ 3 



" For "pertinent" conditions, £ = 500 V/cm, r = 0.5/loti; See 
Eqs (20) and (21). 
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B. Pairwise concurrence of eigenstates 



A quantitative measure of entanglement is provided by the pairwise concurrence function, C\ 2 , 
which becomes unity when entanglement is maximal and zero when it is entirely lacking. The 
general prescription for evaluating C l2 involves somewhat arcane manipulations of the density 



matrix 113911 . However, it becomes simple here as the entanglement arises entirely from off-diagonal 



terms in the Vd-d matrix of Eq.(19). These terms are small, since they are all proportional to y, 
which is < 10~ 5 . Otherwise the off-diagonal terms contain either C x , or C\, factors essentially 
independent of x or i'; for type I qubits, C x ~ 0.4 and for type II qubits C x < 10~ 3 . Accordingly, 
as seen in Tables II and III, the ground eigenstate, *Fi, and the highest excited eigenstate, ¥4, are 
almost solely composed of the basis qubits |00) and |11), respectively, especially for type II. In 
terms of the coefficients in Eq.(22), in this case C i2 is to good approximation just 2d\ or la^, for 
eigenstates 1 and 4, respectively. Thus, for eigenstates 1 and 4, we find 

C 12 = K(x,x'm a /B] (23) 

with weak dependence on x, given by 

K(x) = 0.03752 + 0.003 1 2x + 0.00029* 2 (24) 
and the dependence on jc' is well represented by K(x,x') = \K(x)K(x')~\ l/2 when Ax = x'-x < 10~ 2 . 



The concurrence for a pair of polar diatomic molecules fl29H has this same form (for small Q a /B), 
but the second-order Stark effect makes the K(x) coefficient much larger (> 0.12 for x < 1). 

The C12 function becomes more interesting for the middle eigenstates, y V 2 and TV As seen in 
Tables II and III, for the conditions we refer to as "pertinent" these eigenstates are essentially just 
the |01) and |10) basis qubits, respectively. However, if Ax — » 0, the eigenenergies E 2 and E 3 
become the same. In that limit, even very small y can produce strong entanglement of the |01) and 
1 10) qubits. Figure 2 illustrates how C\ 2 varies as Ax is scanned over a range from well below to 
well above y; at least in principle that can be done by adjusting the 8-field and/or the spacing of 
the dipoles. The curve shown is given by 

C n = 2\a ± \/(l+a 2 ± ) (25) 
11 



with 

(E 3 -E 2 )±\(E 3 -E 2 ) 2 + 4A 2 ] m 
a ± = L_ — (26) 

where A = C|Qq,. This formula for C\ 2 results from omitting all off-diagonal terms in the Vd-d ma_ 
trix except the pair that couple |01) and |10) along the antidiagonal. The eigenstates then become 
*F 2 = ¥ + and ¥3 = V F_, with 

*, = "° > -^ 1> (27) 
Vl + al 

In the limit E 3 - E 2 <sc A (i.e, Ax «c y), where a+ —* +1 and C\ 2 — » 1, the eigenfunctions 
become maximally entangled states, termed Bell states. Figure 2 also displays points obtained 
from numerical diagonalization of the Hamiltonian with all elements included in the Vd-d matrix. 
For both type I (green points) and type II (red points), the numerical results agree very closely with 
the formula given in Eq.(25). It is a striking demonstration of the extent to which matrix elements 
that connect almost degenerate levels generate entanglement. 



C. Inducing large entanglement via resonant pulses 



Under the ultracold conditions needed to localize trapped molecules in the qubit sites, the 
two-dipole system is in its ground eigenstate, ^ ~ |00), wherein the entanglement is very small. 
However, the large entanglement often needed for quantum computing can be induced dynamically 



via resonant pulses to higher eigenstates 11401. 



accomplishing this to use polar molecules in operating quantum logic gates H 



41H . Several procedures have bee n presen ted for 



17 



23 



24 



26 



27 



42144-611 . Here we consider just a rudimentary version, exemplified with the CNOT gate, since our 
chief aim is to compare and contrast the symmetric top qubits of types I and II with the diatomic 



f. fl, 



case treated in ref. Il29ll . 

Figures 3 and 4 give schematic diagrams, analogous to Fig. 10 of ref. [|29|], depicting available 
transitions among the two-dipole eigenstates. Table V lists the corresponding transition frequen- 
cies. In contrast to type I, for type II qubits the contributions from both the rotational constants 
and quadrupole coupling cancel out, hence the transition frequencies depend only on the Stark 
energy shifts and dipole-dipole interaction. Since the entanglement is so feeble for the eigenstates, 
as seen in Tables II and III, for a heuristic description we may speak as if the transitions simply 
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occur between the unperturbed basis qubits. A typical procedure applies a n/2 pulse resonant with 
the transition frequency a>\ to transfer population from the ground eigenstate |00) to the excited 
state |01), thereby putting the system into the state 2~ 1/2 (|00) + |01)). Then a n pulse resonant with 
the transition co 2 between |01) and |11) will put the system into the state 2~ 1/2 (|00) + |11)), which 
is a completely entangled Bell state. The same process can be done applying a n/2 pulse to ct> 3 , 
followed by a n pulse to a) 4 . 

TABLE V: Transition frequencies between eigenstates of two dipoles". 



Type I qubits Type II qubits 
o)i/B 4-y/3-y/6-z/140 x-y/2 
co 2 /B 4-je/3-y/18-z/140 xf + y/2 

cj 3 /B 4-x/3-y/6-z/140 x'-y/2 

C04/B 4-x73-j/18-z/140 x + y/2 

Aco/B y/9 y 



"Here x = fi£/B, y = Q.JB, z = eqQ/B. 

To carry out such procedures, the transition frequencies need to be unambiguously resolved 
from each other. As evident in Table V, for both type I and II qubits, o)\ can be resolved from 
a>2 and 0J3 from 0)4. simply by adjusting the difference in external field strengths, Ax = x' - x. In 
frequency units, a Stark shift of Ax = 10~ 3 for CH3CN is 3 MHz for type I qubits and 9 MHz 
for type II. The relative difference is far more in favor of type II, because oj\ = 35, 869 MHz for 
type I whereas it is only 988 MHz for type II. However, for either type such differences are easily 
resolvable in conventional microwave and radiofrequency spectroscopy. 

Resolving a>i from a>4 and oj 2 from a> 3 presents an experimental challenge. The frequency 
difference is governed simply by the dipole-dipole interaction, since 

Aw = 0)4 - oji = CO2 - CO3 = Q a (.Ci - Cq){C\ - Cq) (28) 

The Aco shift is the essential feature of a CNOT gate: co 3 transfers the target qubit on dipole 1 
from |0) to |1) when the control qubit on dipole 2 is in |0), whereas cl> 2 = 0J3 - Aw transfers the 
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target from |0) to 1 1 > when the control is in |1). For a>\ and co 4 the roles of target and control sites 



are exchanged. Unlike the diatomic case Q29Q , for symmetric tops the cosine elements are nearly 
independent of the external field in the pertinent regime, except via the minor quadrupole terms 
included in Table I. Thus, 

Aoj = Cl a /9 for type I and Aoj = £l a for type II (29) 

Here the significant advantage of type II occurs because both C and C\ are large and of opposite 
sign. In frequency units, for CH 3 CN the Aco shift is only 2 kHz for type I and 18 kHz for type 
EL Again, the relative difference greatly favors type II, since Aco/coi is more than a hundredfold 
larger than for type I. 



As compared with candidate polar diatomic molecules [29J], we expect prospects for resolving 
Aoj for symmetric tops are improved in two ways: (1) The first-order Stark effect enables use of 
a much less strong external field. That should reduce line broadening caused by nonuniformity 
and fringing of the electric field. (2) The choice of Stark components for type II qubits lowers 
the transition frequencies between qubit states down to the radiofrequency range (often factors of 
30-50 lower than transitions between rotational states, which occur in the microwave range). In 
molecular beam spectra, collision free but without trapping in an optical lattice, line widths are 
typically much smaller in the rf region; e.g., 2 kHz or below for Aj = 0, AMj = +1 transitions 



[|36fl . The effect of the optical lattice on line widths is uncertain. It may introduce broadening 



via motional shifts, which are strongly dependent on the well depths required for trapping ||47|1 . 



Such shifts have been avoided for ultracold atoms by use of "magic" optical trapping conditions 



[48], but there might be less scope to do that for molecules. As yet, no line width data have been 
reported for ultracold molecules trapped in an optical lattice and subject to a sizable electric field. 
Thus, although less problematic for type II symmetric top qubits, the feasibility of resolving the 
Aoj shift remains an open question. 

D. Comparison with NMR 

A motivation for considering symmetric top type II qubits is the resemblance t o sp in- 1/2 NMR, 



which has been extensively analyzed in the context of quantum computation |10 

14 



3 SD 

m 
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52J. The 



resemblance stems from the unorthodox choice of +Mj Stark components for type II qubits. That 
renders the effective qubit dipole moments, fi e ff = /j,(cos6), which are essentially independent of 
the external field, equal in size but opposite in spatial orientation. There are further similarities. 
For the generic N = 2 case, the corresponding Hamiltonian for NMR resembles our Eqs. (18) 
plus (19), except for omission of the rotational energy. The molecular dipoles are replaced by 
nuclear spins, the Stark field by a Zeeman field, and the dipole-dipole interaction by spin-spin 
coupling. Thereby our Q, a is replaced by J l2 , the spin-spin coupling parameter. Since the Zeeman 
energy terms are much larger than the spin-spin coupling, the equivalent of our Va-d matrix is 
usually approximated as simply diagonal [3(3]. Accordingly, the eigenstates are then just the basis 
qubits {|00), |01), |10), 111)}, so entirely lack entanglement. That resembles our type II qubits when 
C x = 0, in the absence of quadrupole coupling. 

Another, different sort of similarity arises from the choice of NMR qubits as nuclear spins on 



different atoms within a molecule Q49[|. Even for atoms of the same kind, chemical shifts cause 
the effective external magnetic field to differ at different sites. This corresponds to the role of the 
gradient in electric field, emphasized in Sec.III, wherein Ax > is important both for addressing 
sites and for resolving the |01) and |10) qubit pairs. 

Many procedures for producing dynamical entanglement in NMR systems by means of se- 
quences of radiofrequency pulses have been developed and demonstrated in performing quantum 



gates and algorithms ||10L 



30, 
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■520. The prospects for adapting some of these to polar symmetric 
tops invite systematic study. We will not pursue that here, but mention an example pertinent to 
resolving Aw, the key frequency shift for implementing the CNOT gate. For NMR the analog of 
our Eq.(28) holds, with Aw = 7 12 . 

Even if Aco is too small to be well resolved, another general way to perform a CNOT gate 
has been demonstrated in a NMR spin system Q52D . Because qubits in both sites 1 and 2 are in 
superposition states of |0) and |1), the qubit at site 1 comprises two populations, one coupled to 
the qubit at site 2 in the |0) state and the other to the |1) state there. By means of a n/2 pulse, 
the qubit at site 1 can be rotated into the transverse plane, where both populations will undergo 
Larmor precession, but with different frequencies. After a time ~ 1 /Aoj, the two populations are 
180° out of phase. Then another n/2 pulse can be performed to place both populations at site 1 
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along the z-axis. The net effect is to complete a CNOT gate with the qubit at site 2 controlling 
that at site 1. At least in principle, such procedures, well developed in NMR, seem applicable to 
symmetric top type II qubit states. 



IV. CONCLUSIONS AND PROSPECTS 



The seminal proposal by DeMille fl 1 311 envisioned a quantum computer using as qubits rota- 
tional states of ultracold polar molecules, trapped in an optical lattice, partially oriented in an 
external electric field and coupled by dipole-dipole interactions. Many aspects and variants have 



been extensively studied in the decade since, all considering diatomic molecules H22H29H. As 
the external field has an essential role, the fact that the Stark effect is second-order for diatomic 
molecules has major consequences. The field strength must be sufficiently high to induce exten- 
sive hybridization of rotational states, so that the molecules undergo pendular oscillations about 
the field direction; otherwise rotational tumbling averages out the effective dipole moments in the 
laboratory frame. As discussed in Sec. IIIC, and more fully in ref. I29D, line broadening by the 
high field handicaps resolution of Aa>, the key frequency shift for 2-qubit operations. 

We find that polar symmetric top molecules offer significant advantages. These come primarily 
from the first-order Stark effect, available for all states with K and Mj nonzero. As symmetric 
tops in those states precess rather than tumble, the effective dipole moments are independent of 
the electric field strength (except at high fields). Because there is no need to induce pendular 
hybridization, a considerably lower external field can be used, thereby improving prospects for 
resolving the Aco shift. Moreover, in the first-order Stark effect the ±Mj components are readily 
resolved (not possible for second-order). This enabled considering the \J = l,K = l,Mj = ±1) 
Stark components as the basis qubits (our type II), rather than rotational states (type I). That lowers 
the transition frequencies between eigenstates (cf. Table V) to the radiofrequency range, again 
more congenial for resolving the Aw shift. Even more welcome, the use of ±Mj components as 
qubits brings forth direct correspondences with spin- 1/2 NMR systems. This opens up the prospect 
of exploiting with symmetric tops a wide repertoire of radiofrequency NMR techniques developed 
for quantum information processing. 
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Another prospect for dealing with the small size of the Aw shift involves spatial rather than 
frequency resolution. This is exemplified by quantum computer designs employing supercon- 
ducting flux qubits Il53ll . For these, the generic N = 2 Hamiltonian in the case of transversely 
coupled qubits is much like our Eqs.(18) and (19). Instead of the Stark terms, /i£ and //£', there 
appear single-qubit energy splittings, denoted Ai and A 2 , respectively, and in place of Cl a there ap- 
pears the qubit-qubit coupling energy, denoted by 3 (unrelated to rotational angular momentum or 
NMR spin- spin). The analog of our Vd-d matrix has nonzero elements only along the anti-diagonal 
(equivalent to setting our Co and C\ = 0). However, for typical conditions, 3 « (Ai - A 2 ), the 
analog reduces just to the simple case described under our Eq.(26) and Fig. 2; the correspondence 
replaces our Cj c Q./(E^ - E 2 ) by 3/(Ai - A 2 ). The transitions involved in the CNOT gate (cf. Fig. 
lb of ref. |p3|l ) then occur in degenerate pairs, a>i = a>4 and a>2 = 0)3. Therefore, Aoj = 0, so 
frequency-selective operations are impossible. Yet, one transition of each degenerate pair can be 
selectively suppressed while coherently exciting the other, "by simultaneously driving both qubits 



with the resonant frequency of that pair, employing different amplitudes and phases" 115 311 . This 
method requires spatial resolution sufficient to enable qubits on different sites to be driven individ- 
ually. That may not be feasible for our conditions, with polar molecules separated by only 0.5 /um. 
Such a method is well suited to a proposed design with molecules trapped in QED cavities spaced 
~ 1 cm (!) apart along a superconducting transmission line resonator jlgl. 



As in our previous study of entanglement of polar diatomic molecules Q29U , we provide a 
generic formulation in terms of reduced variables (x, Ax,y,z, w). This makes our results applica- 
ble to a broad class of symmetric top molecules and range of conditions envisioned for proposed 



quantum computers. We also present specific results for the CH 3 CN molecule |38j], regarded as 
a particularly suitable candidate, particularly for type II qubits. Its large dipole moment enhances 
the dipole-dipole interaction and hence the Aoj shift, and its nitrogen atom supplies a quadrupole 
moment that makes the transition dipole Cx nonzero, thereby enabling AMj = +2 transitions 
between the type II qubits. 

Many aspects important for quantum computing with polar molecules are not discussed here 
(trapping operations, sources of decoherence, and much more) because extensive analysis given 



for diatomic molecules 1114 



IMS, 
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2911 pertains as well to symmetric tops. We note an ironic 
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exception. Auxiliary st orag e qubits are sometimes desired to minimize decoherence or to remove 
unwanted information 113011 . Also, "switchable dipole" schemes have been devised to in effect 
turn dipole-dipole coupling "on" or "off' by transferring qubits between states with very different 
dipole moments. For diatomic molecules, such maneuvers typically involve excited electronic 



states; a prototype proposal II 1711 uses CO, for which the dipole moment in the ground X 1 £ + state 
is only 0.1 D, but in the metastable excited a 3 fl state is 1.5 D. For a symmetric top, such things 
can be accomplished more simply by transfer to states with K or Mj zero, where the first-order 
Stark effect vanishes. For example, in the J = 1, K = 1 states of CH^CN under the conditions 
of Table IV, for Mj = the second-order Stark effect [34] yields an effective dipole moment 
of only 0.084 D, whereas for Mj = +1 the first order Stark effect gives an effective moment of 
1.96 D. A transfer Mj = +1 — » 0, without change in the electric field strength, would reduce the 
dipole-dipole coupling 500-fold. 

Symmetric tops offer many other options for qubits. Some, such as hyperfine structure, are 



also available with diatomic molecules. Others are not, such as doublet structures 115411 produced 
by tunneling through barriers to inversion (e.g., in NHt,) or internal rotation (e.g., in CH3CF3). If 
inversion is fast (~ 1 Hz for NH 3 in ground state), the dipole flips rapidly and the Stark effect is 
second-order, whereas if inversion is slow (e.g. ~ 1 year for AsH{), it is first-order. For internal 
rotation involving a three-fold barrier, the tunneling doublets occur as a nondegenerate A state, 
and a doubly degenerate E state; the Stark effect for A is second- order, for E first-order. 

For both diatomic and symmetric top molecules, under conditions considered amenable for 
proposed quantum computers, the entanglement of eigenstates and the associated pairwise concur- 
rences are very small. Furthermore, it is not needed in the eigenstates, because the entanglement 
required for computations is actually induced dynamically. The role of dipole-dipole coupling as 
the source of eigenstate entanglement, via the off-diagonal terms of Eq.(19), therefore is irrelevant. 
Its important role is determining a different eigenstate property, the Aco shift, via Eq.(28). The eval- 
uation of Aw does not require eigenfunctions, only eigenvalues. This is a liberating perspective in 
considering analysis of multidipole systems well beyond N = 2. 



Mindful of the somewhat metaphysical status often accorded to entanglement B55Q , we men 



tion that fundamental theory shows that even for symmetric tops, the "true molecular eigenstates 
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should not have first-order Stark effects" Q56Q . That is because the full permutation-inversion group 
for a molecule shows that the only levels allowed by quantum statistics are nondegenerate. Yet 
both theory and experiment confirm that a quasi-first-order Stark effect does appear in the pres- 
ence of even a very weak field (< 0.3 V/cm) that introduces coupling between nearly degenerate 
states. Hence, the very existence of first-order Stark effect in molecules comes from field-induced 
entanglement. 
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APPENDIX A: QUADRUPOLE COUPLING 



As outlined in Sec.IIB, we use the "strong-field" approximation 11321. 13711. appropriate when the 
Stark shifts are much larger than hyperfine splittings introduced by quadrupole coupling. We need 
to evaluate contributions from H e to be added to the qubit eigenvalues of Eqs.(10) and (13). Also, 
we need to obtain, by diagonalizing H 5 + H e , the modified qubit eigenfunctions that arise from 
mixing of the Mj Stark components with the Mj nuclear spin components. These are required 
to determine the quadrupole contributions to the cosine elements of Table I. The requisite matrix 
elements of the H e Hamiltonian, 



<7, K, I, Mj, Mj\H q \J, K, I, M'j, M'j) 



(Al) 



are given in Eq.(33) of ref. |37|]. All contain a common factor, 

eqQ / 3K 2 



P(J,K,I) = 



4(27 - 1)(27 + 3)(2/ - 1) \ 7(7 + 1) 



(A2) 
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For the qubit states we consider, 



P(l, 1, 1) = eqQ/40 and P(2, 1, 1) = -eqQ/m 



(A3) 



The elements of H e comprise a 9 x 9 matrix labeled with Mj = 1, 0, -1 and Mj = 1, 0, -1. The 
first order energy of the quadrupole hyperfine components is given by the diagonal elements, 



E Q = P(J, K, I) [3Mj - J (J + 1)] [3Mj - 1(1 + 1)] 



(A4) 



Because the sum M F = Mj + M/ is a good quantum number, the matrix is block diagonal, with 
five submatrices corresponding to M F = 2, 1,0, -1, -2 (respectively 1 x 1, 2 x 2, 3 x 3, 2 x 2, 
lxl). We deal only with the M F = block, containing elements connecting the (Mj,Mi) = 
+ 1, -1; 0, 0; and - 1,+1 hyperfine components: 



P(l,l,l) 



( \ 
1 -3 6 



-3 4 -3 
6 -3 1 



(A5) 



-1 - V3 6 
-V3 4 -V3 
6 -V3 -1 

To the diagonal elements of these matrices, we add the Stark components, from Eq.(6), E s = 
-QuE/J(J + l))Mj, then carry out diagonalization to obtain the M F = eigenfunctions. As speci- 
fying Mj automatically specifies Mj, we denote the eigenfunctions simply by ^(7, Mj), expressed 
as linear combinations of the basis functions (f>(J,Mj). Here we revert to wavefunction notation, 
to avoid confusion with the bra notation used for qubits. Also in labeling the eigenfunctions, 
we adorn Mj with a tilde, to indicate it is no longer a good quantum numbers because the Stark 
and spin states are mixed. Performing numerical diagonalizations led to recognition that, for 
/j£» eqQ, the eigenfunctions are well approximated using for each J a single mixing coefficient; 
for 7 = 1: 

¥(1,-1) « (1 -a 2 )(f>(l,-l)- 00(1,0) + a0(l,+l) (A7) 



¥(1,0) *a0(l,-l) + (l -a 2 )0(l,O)-a0(l,+l) 



(A8) 



20 



¥(!,+!) * -a0(l,-l) + a0(l,O) + (1 -a 2 )0(l,+l) 



(A9) 



and for 7 = 2: 



¥(2, -1) * (1 - 2£ 2 )0(2, -1) + b<f>(2, 0) - V3£0(2, +1) 

¥(2,0) w -^0(2, -1) + (1 -Z7 2 )0(2,O) + ^(2,+l) 
*F(2,+i) « V3&0(2, -1) - ^0(2, 0) + (1 - 2Z? 2 )0(2, +1) 



(A10) 
(All) 
(A12) 



The coefficients a and b are small positive numbers, determined by w = \eqQ\//uE. From our 
numerical results, we find 



a = 0.1522w and b = 0.1789w 



(A13) 



These values are accurate within 1% for w < 0.1. Since Mj = +1 for our qubit states, as defined in 
Eqs.(9) and (12), we now specify them further as the hyperfine components *¥(J, -1) and ¥(7, +1); 
thus for type I, 



|0> = \J= l,Mj = -\,Mi = +1) and |1) = \J = 2, Mj = -1, M l = +1) 



(A14) 



and for type II, 



|0> = |7=l,Af 7 = +l,Af/ = -l> and |1> = \J = 1, Mj = -1, M l = +1) 



(A15) 



The quadrupole terms in the cosine elements of Table I result from using the mixing coefficients 
of Eq.(A13) with Eqs.(A7) and (A10) for type I qubits and Eqs.(A7) and (A9) for type II together 
with Eq.(7) of Sec.II. In particular, for type II this gives C x ~ a{\ - a 2 ). 

Symmetric top molecules, other than CH 3 CN, which contain one nucleus with spin / = 1 on 
the symmetr y ax is, include: NHi, and NF$, where 14 Af has eqQ = -4.09 MHz and 7.07 MHz, 



respective 



tively |36, 



y|57j 



5811 : and C// 3 D and CF 3 D, where 2 D has eqQ = 191 kHz and 171 kHz, respec- 



59J]. In many halide molecules, such as CH^X, the halogen nuclei have / > 1 and large 



quadrupole coupling constants 



weak-field approximation Q37 



60] 



Treatment of such cases requires use of an intermediate or 
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APPENDIX B: ENTANGLEMENT OF TWO DIPOLES 



For the ranges of reduced variables specified in Sec. Ill: x < 1; 10~ 4 < Ax < 10" 2 ; y < 10" 5 ; 
\z\ < 5 x 10~ 3 ; w < 0.1, we have obtained explicit formulas for the coefficients of the basis qubits in 
Eq.(22), {a,-, bj, q, d { }, that determine the two-dipole eigenstate entanglements. Tables VI and VII 
give these formulas for types I and II qubits, respectively. Also included are corresponding values 
of the pairwise concurrence, C i2 , for the eigenstates; these conform well to the approximations 
of Eqs.(25) and (27). The corresponding eigenvalues and orders-of-magnitide of the coefficients, 
under conditions listed in Table IV, are in Tables II and III. Contributions from quadrupole coupling 
are not included in Table VI because these only slightly affect the entanglement for type I qubits. 
The quadrupole contributions are included in Table VII because for type II qubits these are the sole 
source of eigenvalue entanglement (since without them C x = and the Vd-d matrix of Eq.(19) is 
diagonal). The quadrupole contributions enter the entanglement coefficients in various powers of 
the ratio of the quadrupole coupling to the Stark energy, w", ranging from n = 2 to 4. In the 
concurrence values, the same dependence on w n appears. 

Tables VI and VII both pertain to the regime Ax » y, where the Stark shift between adjacent 
qubit sites is much larger than the dipole-dipole interaction. At present, this regime appears most 
relevant for implementation prospects. As illustrated in Fig. 2 and Eq.(27), therein the eigen- 
functions differ little from the basis qubit states, and entanglement is slight. The extreme opposite 



limit, Ax = 0, has been analyzed in ref. [|23D; there the eigenfunctions X P 2 and ¥3 become the 
maximally entangled Bell states, 2 _1/2 (|01) ± |10)). An interesting consequence emerged. For 
operation of the CNOT gate, it was concluded that a preliminary pulse of bandwidth much wider 
than the dipole-dipole interaction should be applied. It would entirely undo the entanglement by 
forming + combinations of the Bell states and thereby unwed the |01) and |10) qubits. 
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TABLE VI: Eigenfunction entanglement coefficients, type I qubits' 



*F; = a;|00> + 6,101 > + cj|10> + d/|ll) 
fli - (0.07.y) 2 /2; h = c x = (0.048 + 0.0044x);y; 

d x = -04 - (0.019 + 0.0017jc)y; C i2 - ^(x, x').y « 2di ; 
a 2 = -(0.028 + 0.0026x).y; 62 = 1 - c 2 /2; c 2 = -0.454);/ Ax; 

J 2 - (0.009 + 0.0009jc)y; C 12 « 2|c 2 |; 
a 3 = -(0.062 + 0.0058x);y; fc 3 - -c 2 - 0.454;y/Ax; c 3 - b 2 \ 

d 3 = (0.021 + 0.0019x).y; C 12 « 2|c 2 |; 
a 4 = -rfj = (0.019 + 0.0017x)j; b 4 = c 4 = -(0.016 + 0.0015x)j; 
d 4 = -a A = 1 - (0.04;y) 2 /2; C n = K(x, x')y * 2a 4 ; 

"Coefficients {a,, fcj, c,, c/,} of eigenfunctions ; = 1 — > 4 obtained from numerical diagonalization of the matrices 
of Eqs(18) and (19). Table II gives the corresponding eigenvalues as well as orders-of-magnitude of the coef- 
ficients under conditions listed in table II. Values are included for the pairwise concurrence, C\2, and conform 
well to the approximations of Eqs.(25) and (27). 



TABLE VII: Eigenfunction entanglement coefficients, type II qubits" 

l F ; = a i |00> + ^|01> + c,|10> + J ; |ll> 
a\ = 1 - b\\ bi = d = -0.0786wVz; 

d\ = -0.01 18wVz; C12 = 0.0247w 3 y/z * ld\ ; 

a 2 = d 2 = a 3 =d 3 = 0.0743w 2 j/z; b 2 = -c 3 = -0.0225w 2 y/Ax; 

c 2 = 1 - a\\ d 2 = a 2 ; C\ 2 = 0.044w 2 y/Ax « 2|& 2 |; 
a 3 = d 3 = a 2 = d 2 \ b 3 = c 2 = 1 — a^', c 3 = —b 2 - 0.0225w 2 y/Ax; 

d 3 =a 3 ; C X2 = 0.044w 2 y/Ax « 2|c 3 |; 

a 4 = 0.0107wVz; b 4 = c 4 = -0.0715w 2 y/z; 

d A = 1 - b\\ Cu = 0.0205 w 3 y/z * 2a 4 ; 

"Footnote to Table IV pertains have as well, except that corresponding eigenvalues and order-of-magnitude 
values are in Table III. Contributions from quadrupole couplings are included with w = \z\/x = \eqQ\/(iE and 
z = eqQ/B 
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FIG. 1: (Color online)Stark states for a polar symmetric top molecule, as functions of /J.E/B. (a) Eigenen- 
ergies for Mj = and ± 1 components of K = 1 levels for / = 1 and 7 = 2 states and (b) corresponding 
expectation values that determine effective dipole moments, fx e ff = (cos). States used as basis qubits are 
labeled |0> and |1): type I (green) are Mj = -1 for J = 2 and type II (red) are Mj = +1 and -1 for J = 1. By 
virtue of the ordinate scale used, (a) as well as (b) applies to any symmetric top molecule (treated as rigid, 
without fine or hyperfine structure). 
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FIG. 2: (Color online) Pairwise concurrence C\i for eigenstates 2 and 3 of two symmetric top dipoles 
entangled via dipole-dipole interaction, as a function of the ratio of the difference of the eigenvalues, (£3 - 
E%), to the element, C\ that connects the |01) and |10) basis qubits in the Vd-d matrix of Eq.(19). The 
difference (£3 - E-i)lB = Ax/3 and Ax, for type I and type II qubits, respectively, as seen in Tables II and 
III. Points (green for type I, red for II) were obtained from numerical calculations including all elements 
of the Vd-d matrix; curve (black) from the minimalist 2x2 model of Eqs.(25-27). The same C\% function 
applies to spin-1/2 NMR systems, with £3 - £2 = g^N(H' - H) and C|f2 a replaced by 5/12, where g is 
the nuclear g-factor, [i^ the nuclear magneton, H the magnetic field strength, and Ju the spin-spin coupling 
constant. 
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FIG. 3: (Color online) Schematic energy levels for type I qubit eigenstates of two symmetric top dipoles. At 
left are indicated qubit basis states, with corresponding eigenenergies from Eqs.(18) and (19). Contributions 
from quadrupole coupling are not shown (but included in Tables II and V). At right are transitions that are 
involved in CNOT operation: a>i transfers the dipole at site 2 from |0) to |1), with dipole at site 1 remaining 
in |0), then 2 transfers dipole at site 1 from |0) to |1) with dipole at site 2 remaining in 1 1 > . The same result 
could be reached by co^ followed by 0)4. Transition frequencies are given in Table V. 
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FIG. 4: (Color online) Schematic energy levels for type II qubit eigenstates of two symmetric top dipoles; 
format as in Fig. 3. Eigenenergies, including quadrupole coupling are (not shown) are given in Tables III 
and transition frequencies in Table V. 
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